Disordered d-wave superconductors with chiral symmetry 
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A two-dimensional lattice model for d-wave superconductor with chiral symmetry is studied. 
The field theory at the band center is shown to be in the universality class of U(2n)/0(2n) and 
U(2n) nonlinear sigma model for the system with broken and unbroken time-reversal symmetry, 
respectively. Vanishing of the beta function implies extended states at the band center. Density of 
state vanishes as a cubic function of the energy at the band center for the former case, while linear 
for the latter. 
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Symmetries play a crucial role in random critical phe- 
nomena. The rotational and time-reversal invariances 
specify the well-known universality classes whose under- 
lying symmetries are orthogonal, symplectic, and uni- 
tary ||. The variety of universality classes, however, has 
turned out to be richer. Above all, Altland and Zirnbauer 
have recently discussed universality classes for quasi- 
particles in disordered spin-singlet BCS Hamiltonian, de- 
noted in their article by C and CI for the case with bro- 
ken and unbroken time-reversal symmetry, respectively. 
These universality classes are actually involved with vor- 



L quasi-particle trans- 
| , and the spin quan- 



tices in s-wave superconductors 
port in d-wave superconductors I 

turn Hall (QH) transitions |7). Especially, dirty d-wave 
superconductors have attracted a lot of interest, because 
they can be described by Dirac Fermions || due to gap- 
less quasi-particle spectrum. Random Dirac Fermions 
actually give nontrivial critical points |9) . The criticality 
of the density of state (DOS) of disordered d-wave super- 
conductors is, for example, still controversial |j 10 f|. 

Another possibility is the class for the two-sublattice 
model pd| ], where Hamiltonian has a special symmetry 
which will be specified momentarily. This symmetry will 
be referred to as chiral symmetry |l2| . One example in 
two-dimension is the random flux model [fl3f and another 
is the random hopping fermions with 7r-flux [fl4|| , both 
of which have an isolated delocalized states at the band 
center. Various universality classes, including those men- 
tioned above, are well summarized by Zirnbauer 

In this paper, we study a system with chiral symme- 
try in d-wave superconductors. The basic idea is that 
the pure d-wave lattice system is a kind of two-sublattice 
model. It is, therefore, natural to include randomness 
keeping the same symmetry of the pure Hamiltonian. 
The model we study has a SU(2) symmetry as well as 
the chiral symmetry. It is also interesting to study a 
system with such enhanced symmetries. We show that 
the field theory of the model with broken and unbro- 
ken time-reversal symmetry is in the universality class of 
U(2n)/0(2ra) and U(2n) nonlinear sigma model (NLSM), 
respectively. It is interesting to compare the classes to 
those of two-sublattice models without SU(2) symme- 
try, which are in the class of U(n) and U(2n)/Sp(n) 
NLSM or to those of d-wave superconductors without 



chiral symmetry, which are in the class of Sp(n)/U(2n) 
and Sp(n) NLSM for the case with broken and unbro- 
ken time-reversal symmetry, respectively. We conclude 
that regardless of the time-reversal symmetry the present 
model has extended states at the band-center, which are 
associated with diffusive spin transport Hd] ■ We discuss 
the behavior of the quasi-particle density of state (DOS), 
which should vanish at the band center as E 3 and E in 
broken and unbroken time-reversal cases, respectively. 

Let us start from the lattice Hamiltonian of a d-wave 
superconductor defined on the square lattice in two- 
dimension lyi, 



H 



£4 



(1) 



Hermiticity and SU(2) symmetry (spin-rotation) require 



Uj = and Aij — Aj i: respectively. In the absence 



of randomness, the pure d-wave Hamiltonian is given by 
tj,j±x = tj,j±y — ~to, ^j,j±£ — —Ajj±y — Ao, and oth- 
ers are zero, where x = (1,0) and y = (0, 1). It is easy 
to see that the pure Hamiltonian change the sign under 



the transformation 



(-) 



which is the chi- 



ral symmetry as mentioned-above. Therefore, the pure 
system has the SU(2), time- reversal, and chiral symme- 
try. By introducing disorder which breaks some sym- 
metries explicitly but keeps the others, the model shows 
various kinds of universality classes. Those studied so 
far are as follows: If we keep only SU(2) symmetry, the 
(replicated) model has Sp(n) symmetry, which is spon- 
taneously broken to U(2n). If we keep in addition the 
time-reversal invariance, the model has an enhanced sym- 
metry Sp(n)xSp(n), which is broken to Sp(?i) (^](|. In 
this paper, the symbol (i, j) in Eq. (Q) is restricted to the 
nearest neighbor pairs, which keeps the chiral symmetry, 
and universality classes of the d-wave superconductors 
unknown so far are investigated. 

The pure Hamiltonian has four nodes, where gapless 
quasi-particle excitations exist p). To study the low- 
energy properties of the system governed by them, let us 
take the continuum limit near the nodes B, 
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Cji/a~i i ' +i »xlu( x )+ i 31 Jv x\ 2 i( x ) 

+i-j*+3v x \ 12 (x) + P*-iy x \ 22 (x), (2) 

where a is a lattice constant, x = aj, and indices a, a, 
and i of the field Xaa i are associated with spin, particle- 
hole, and nodes, respectively. Namely, the field x at each 
x lives in the space V = C 2 ® C 2 ® C 2 . The pure Hamil- 
tonian ||[| is then H p = J d 2 x X ^H p Xi where 



Hv — I2 



(x <-> y) I ' 



(3) 



and the coordinates have been transformed as x, y — > 
±x ^ ■ Here the explicit matrix in Eq. (||) denotes the 
node space, and gamma matrices denotes the particle- 
hole space, defined by 71 = ra^ and 72 = r~ 1 o\ with 
r = y/vp/vA, and v = ^/i>fi>a, where vp ~ 2\[2ta and 
v A = 2 % /2A . 

The symmetries of the lattice model are translated into 
the continuum model as 



cue- 1 = -n\ 
vnv~ l = -n, 
THT- 1 = -n, 



C = ia 2 ® I2 ® I2, 
V = i<r 2 ® ia 2 ® I2, 
T = CP, 



(4) 



which describe, respectively, the SU(2), chiral, and time- 
reversal symmetry. The total Hamiltonian density is 
given by Ji = 7i p + Hd, where Hd is disorder poten- 
tials satisfying Eq. (Q). In what follows, we study the 
model with broken and unbroken time-reversal symme- 
try at the same time. For the former case, the third 
condition should be omitted. It is easy to write down ex- 
plicit disorder potentials using these conditions, though 
it is not necessary in the following calculation. 

To study the DOS and the conductance of the spin 
transport pa], let us define the Green functions, 



G(x) =tr v (x\{ie-H)- 1 \x), 
K(x,x') = try |(x|(ie-H)-V)| 



(5) 



where try is the trace in the V space. Although it 
may be easy to introduce the generating functional of 
these Green functions, we have to prepare some nota- 
tions. Firstly, the generating functional expressed by 
path-integrals over Fermi fields is defined in a standard 
way by introducing replica, Xi -> Xia and X i ~* Xai, 
where i and a are indices denoting V and replica space 
Wr = C™, respectively. Lagrangian density is then C = 
— triy R X — ^0 Xj where triy R is the trace in the replica 
space. Moreover, we introduce an auxiliary space || to 
reflect the symmetries in the V space to an auxiliary field 
introduced later [See Eq. (EM W R W = W R <g> W A , 



with Wa = C 2 ® C 2 ® C z , where each space in Wa is 
associated with the SU(2), chiral, and time-reversal sym- 
metry, respectively, and in this extended space, Fermi 
fields are denoted by ip a i and ipia, which are subject to 



V> =CV' t 7" 1 J 
ip = iPipr -1 , 
ip = iTipir^ 1 . 



(6) 



Matrices 7, r, and it are defined in the W space, given 
by 7 = l n ® i& 2 (8 I2 <8 I2, t = 1„ (8 1 2 (H)^ " 2 ®J-2> and 
7r = 1„ ® I2 ® I2 ® fi- The identity tr M /(oj-0V ; ~ ipHip) — 
tr^y R x(«e — H) X , where w — iel n ® 03 ® 03 03 , leads to 
the generating functional Z = jT)ipT)ipe~ s with 

S = - Ji d 2 xtr w (Wv ~ ■if'Hip + J?AV) • (?) 
Green functions are expressed as 



G(x) = d WA lim ((i/jip) 11 (x)) , 

n— »0 \ / 



K(x,x') = -d 2 Wh lim l ( W)ia(a/)(^)ai(x)) , (8) 

where 1 = a<g>T®T<g)T and 2 = /3(g) J. ® f ® j, and d^y A 
is the dimension of the auxiliary space, given by 4 and 
8 for broken and unbroken time-reversal case, respec- 
tively. Assume that disorder potentials obey the Gaus- 
sian distribution J VHd exp(— J d 2 xtrvH 2 ). Then it 
is easy to average over them, as usual. Here, some 
comments may be useful Q: Firstly, disorder potentials 
are integrated out by using — ^tryH^ + trvH-dipip = 
— j^itryTid — gipijj) 2 + §try(V>V0 2 - It turns out that 

the integration over Tid is automatic because tpip sat- 
isfy the same symmetries as those of Hd due to Eq. (|^). 
Secondly, try (ipip) 2 = — trw (ipip) 2 , which is converted 
into Yukawa-type interactions by introducing an auxil- 
iary matrix field defined in the W space into the action 
in the form j^tr w(Q + gtyty — 0S) 2 . Integrating out the 
Fermi fields, we end up with an effective action, 



S= -^-Jd 2 xtr w (Q 2 -2Qu 



InDetyg.iy (1 <g) Q - H p <8> 1) 



(9) 



Here we have set J = for simplicity and the anti- 
Hermitian auxiliary field Q = —Q' is subject to 

Q = -7QV\ Q = -tQt-\ Q^-nQn- 1 . (10) 

To get a saddle-point solution, set w = 0, and assume 
Q is diagonal and spatially constant. Then, saddle-point 
equation reduces to 



g d WA 



try(Qo-Hp) 1 (x,x), 



(11) 



which gives Qo = ivo^-n a 3 ® &3 ® U3 as a solution. Here 
vq ~ A exp( — irdw A v 2 / Ag) with A being high-energy cut- 
off. This solution gives rise to an exponentially small 
density of state at the band-center. 
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The effective action has U(2n) and U(2n)xU(2n) sym- 
metry for the model with broken and unbroken time- 
reversal symmetry, respectively. To be concrete, let us 
consider the former case. For the action function to be 
invariant under the group action Q — > gQg -1 , g should 
satisfy 7 = .973* and rgr^ 1 = g. They are explicitly 
given by 

9 = v( U _)vt, Q = v(_ qf (12) 

where matrices in Eq. ( |l2| ) denote the chiral space of W, 
and u is 2n x 2n unitary matrix in the replica and SU(2) 
space of W, while q is a complex matrix with a condi- 
tion q t = q. V is a matrix defined by V = V\V% with 
V\ = l n ®/i®/z and V 2 = \ n ® I2 ® + 03) + In® 02 <S> 
5(12 — C3)) where /1 = -^75(03 +CT2). It is, therefore, easy 
to show that the saddle-point solution is invariant under 
the group action if u l u — 1. Thus, the U(2n) symmetry 
is spontaneously broken to 0(2n), and the saddle-point 
manifold is U(2n)/0(2n). 

Local quantum fluctuation around the saddle-point 
yields Goldstone mode, which governs the low energy 
properties of the present system. To derive an effective 
action of this mode integrating out the massive mode, let 
us decompose the field Q into transverse and longitudinal 
modes as Q(x) = T(x)(Qo + L(x))T> (x), where 



t v( u v ^v\ L = v^_ p p )yt. (13) 



Here U(x) € U(2n) /0(2n) and P(x) is a real matrix with 
P t = P. Next transformation properties of U(x) and 
P{x) fields should be examined. The action of global u £ 
U(2n) induces uU(x) = U'(x)h(u,U(x)), where U'(x) e 
U(2n) and h(u,U) G 0(2n). It should be noted that 
h(u, U) is a nonlinear function of U(x) dependent on x. 
Therefore, we have the following transformation laws of 
U and P, U -> uCZ/i" 1 and P -> frPft" 1 . 

Keeping this in mind, we proceed to calculate the ac- 
tion of the Goldstone mode. The point is the derivative 
expansion of the In Det term in Eq. d9h , 



lnDety®w(l <8> Q - Hp ® 1) 

1 2 

~ Try^^A^G — -Tiv®w (ApG) . 



(14) 



where Ap is the free propagator defined by A^ 1 = 
1 ® Qo - Hp ® 1 and G is 



G= (1 <g> Tt)(W p <g> 1)(1 <g> T) - W p ® 1 
(x «-> y) 



(15) 



The explicit matrix in the above denotes the node space 
in V. It is stressed that the effective action should be in- 
dependent of the gauge associated with the local 0(2n) 



transformation. To see this, let us write down the field 
T^dfjT in W space, 



T^T = V 



V u + A. 



Vu, - A, 



(16) 



The transforma- 



hAph 1 and 



where A^V^ = ^Wd^U T XJ^d^U). 
tion laws of these fields are A^ 

hV^hr 1 + hd^hr 1 , which tells that is a gauge field as- 
sociated with the hidden local 0(2n) symmetry Jl||. In 
the leading order of the derivative expansion, therefore, 
only the gauge covariant A^ field appears in the action. 

By using these, manifestly gauge-invariant action of 
the Goldstone mode can be obtained. The principal term 
is calculated as Tt v&w (A F G) 2 = -|Tr WRS A^. Af- 
ter gauge- fixing C/' = U, we have 



S = I d 2 xtTw R 



v J*(u 2 + u- 2 ) 



(17) 



Here Wrs denotes a part of the space W restricted to the 
replica and SU(2) spaces omitting the chiral space, and 
1 _ 1 "f+"a _ This coupling constant is just the same as 

O 47T VFVA i— i 

that derived by Senthil at al. |Q . So far we have obtained 
an effective action for broken time-reversal system. It is 
easy to follow similar calculations for the unbroken time- 
reversal system. The results are the same as Eq. ( |l7| ) 
including the coupling constant, but with U(x) £ U(2n) 
and a breaking term 2 ^tr H / RS 3?(?7 2 + U~ 2 ). 

In the process of the renormalization, it is shown that 
the operator fTr 2 ^^ = -LTr^t/ 2 ^- 2 is needed. 
In the leading order we have taken above, the coupling 
constant is - = 0. However, higher loop expansion of 
the In Det term actually gives a finite coupling constant. 
After the replica limit n — > 0, the renormalization group 
equations at one-loop order are given by p?| , pT[ | 



db 
dl 



= 0, 



dc 
dl 



-OL\C 



dc , , If b 2 



(18) 



where (ai,/3i) = (|, §) and (±,0) for U(2n)/0(2n) and 
\J(2n) models, respectively. The spin conductance is re- 
lated only with the coupling constant b as r = 2ai7r<7 s , 
which is obtained by the bare diffusion constant calcu- 
lated from Eq. (||). Therefore, vanishing of the beta 
function for b ]l8| implies delocalized state at zero energy, 
which is due to the chiral symmetry. On the other hand, 
the behavior of the DOS depends on the coupling con- 
stant c. Since e has the same dimension of the energy, we 
can compute, according to Gade pJ| , a rough estimate 
of the DOS for finite E as p{E) ~ Ag-V^ib 2 )- 1 in(A/B)_ 
This yields an enhancement of DOS near the zero ener- 
gies. However, the result from perturbative calculations 
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needs alternative consideration for DOS at the zero en- 
ergy. 




FIG. 1. Schematic illustration of the density of state as a 
function of energy. Dashed-line denotes DOS of pure system, 
while dashed-dotted-line (full-line) denotes DOS for broken 
(unbroken) time-reversal symmetry. The bump means the 
enhancement obtained by the renormalization group analy- 



Since derealization occurs at the zero energy, the lo- 
calization length is quite long, i.e., the order of the system 
size L near the zero energy. Accordingly, in the energy 
scale much smaller than the Thouless energy jf, where 
D s is a spin diffusion constant, quasi-particles are diffus- 
ing all around the system keeping the symmetry of the 
Hamiltonian, and hence, the spatial dependence of the 
system is smeared out and we can describe the system 
by a random matrix theory. Taking only zero mode of Q 
into account, we actually have from Eq. (^|) an effective 
action of a random matrix theory with the symmetries 
((§ but defined in V = C N ® C 2 ® C 2 ® C 2 , where 
N = L 2 . The Hamiltonian Tt is now not a field but a 
quantum mechanical one subject to Eq. (Q). To cal- 
culate the DOS near the zero energy of such a random 
matrix ensemble, it may be convenient to rotate the basis 
by a orthogonal transformation and to switch into more 
convenient basis. We can choose C — In <S> I2 ® i&2 ® I2, 
V = In <£> I2 ® I2 <8 03. Let us omit the SU(2) space 
in the above, because it is irrelevant. A little thought 
tells that the Hamiltonian describes the tangent space of 
Sp(W A)/Sp(A0xSp(A0 and \](N, N)/U(N)x\J{N) for 
the case with broken and unbroken time-reversal symme- 
try, respectively, and accordingly, belongs to CII (chiral 
GSE) and Alll (chiral GUE) in Zirnbauer's classifica- 
tion p]. We can now diagonalize the Hamiltonian as 
H = Udi&g(8, —8)U', where U is a unitary matrix. The 
Jacobian of the change of variables from Hij to f/y and 9i 
is J = const.U l 9f 2 n i<j {9 2 ~9 2 )^ 1 where (a 2 ,/3 2 ) = (3, 4) 
and (1, 2) for broken and unbroken time-reversal case, re- 
spectively. /3 2 describes the level repulsion of eigenvalues 
each other, while a.i describes the repulsion of each pair 
9 and —9 and determines the DOS for small 9 M. Based 
on this argument, we expect that DOS vanishes at zero 
energy as E 3 and E, respectively, for the system with 
broken and unbroken time-reversal symmetry. Interpo- 



lating the results, we expect the behavior of DOS as Fig. 
[I]. It is quite interesting to test the present conjectures 
by numerical calculations. Especially a cubic dependence 
of DOS on the energy has not been known so far in dis- 
ordered d-wave superconductors. 
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